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ABSTRACT 

We report on the first high frequency VLBI observations of the nearby broad 
absorption hne quasar (BALQSO), Mrk 231. Three epochs of observations were 
achieved at 15 GHz and 22 GHz, two of these included 43 GHz observations as 
welL The nuclear radio source is resolved as a compact double. The core com- 
ponent experienced a strong flare in which the flux density at 22 GHz increased 
by > 150% (45 mJy) in three months. Theoretical models of the flare imply that 
the emission is likely enhanced by very strong Doppler boosting of a highly rel- 
ativistic ejecta with a kinetic energy flux, Q ~ 3 x iC'^ergs/sec. Combining our 
data with two previous epochs of 15 GHz data, shows marginal evidence for the 
slow advance of the secondary component (located ~ 0.97 pc from the core) over 
a 9.4 year span. We estimate that the long term time averaged kinetic energy 
flux of the secondary at Q ~ 10*^^ ergs/sec. Low frequency VLBA observations 
indicate that the secondary is seen through a shroud of free-free absorbing gas 
with an emission measure of ~ 10^cm~^pc. The steep spectrum secondary com- 
ponent appears to be a compact radio lobe that is associated with a working 
surface between the ram-pressure confined jet, and a dense medium that is likely 
to be the source of the free-free absorption. The properties of the dense gas are 
consistent with the temperatures, displacement from the nucleus and the column 
density of total hydrogen commonly associated with the BAL wind. 
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Introduction 



At a redshift of 0.042, Mrk 231 is one of tlie nearest radio-quiet quasars to earth. 
Mrk 231 is a well studied object because it is one of the brightest AGN in the infrared 
sky. It is a quasar that is partially obscured by dusty gas as evidenced by its red op- 
tical/UV spectrum (ILipari et all 1 19941 : ISmith et all Il995l ). The host galaxy has a strong 
starburst contribution, but the obscured AGN bol ometric luminosity , Lfep;, can be inferred 



from a pronounced 10 fi peak in the IR spectrum (IFarrah et alll2003l : iLonsdale et alll2003l ). 
Using the extracte d IR AGN luminosi t y, we estimate from QSO composite spectra (see 
Punslv and Tingavl J2005h : IXelfer et al I J2002h : Izheng et all Jl997h ) that the total thermal 



luminosity created by the central engine is Lf,oi ^ 2 x lO^^ergs/s and Mrk 231 would appear 
as a quasar with My ^ —25.3 if not f o r the semi-transparent dusty shroud (in agreement 
with the estimate in iBokensberg et al (|l973)fl This equates to about 2.7 magnitudes of 
visual extinction created by th e dusty enve l oping gas. Comparing the composite spectrum 
to the optical/UV spectrum in ISmith et all (Il995l ) indicates even more severe reddening in 
the blue band, 3.1 magnitudes at 4400A. It is the closest radio quiet quasar in which there 
is sufficient radio flux to permit high signal to noise imaging with high frequency VLBI 
(Very Long Baseline Interferometry) observations. As a result of the ability to utilize the 
high resolution of VLBI, this nucleus can be resolved with higher precision than any other 
radio quiet quasar nucleusE Our 43 GHz observations can resolve the nucleus to within 
5 X lO^^cm (~ 0.16 pc). It is a fortuitous circumstance that this object is a broad absorption 
line quasar (BALQSO) as well. About 10% - 15% of quasars show broad absorption line 
(loosely defined as absorbing gas that is blue shifted at least 5,000 km/s relative to the QSO 



rest frame and displaying a spread in velocity of at least 2,000 km/s) outfiows (IWeymann 



19971 ). It is widely believed that most or all radio quiet quasars have BAL outfiows, but the 
designation of BALQSO depends on whether the line of sight intersects the 



^Throughout this paper, we adopt the foUowing cosmological parameters: Hq=7\ km/s/Mpc, 0,^ = 
0.73 and Jl,„ = 0.27 in all our calculations. At z = 0.042, 1 mas angular separation equates to 0.818 pc 
displacement in the sky plane in the adopted cosmology. 

The radio loudness, R, is usua Uy defined as a 5 GHz fl ux density 10 times larger than the 4400A flux 



density, R = '^SGHz/S'^^^^j^ < 10 ( Kellermann et al 19891 ). This definition is not appropriate, if there is 



intrinsic extinction as in BALQSOs. One must correct f or the intrinsic blue band extinction and recompute 
the R with dereddened blue fluxes, ([Becker et a i l2oon[) . Noting the blue band extinction in the text and 
using the historic 5 GHz flux densities from NED yields a dereddened radio loudness, 1.4 < i? < 3.8. This is 
close to the ra dio intermediate q uasar category, 10 < i? < 250, which are often just Doppler boosted radio 
quiet quasars ( Falcke et all 1996 ). We show that Doppler boosting greatly enhances R in Mrk 231 as well. 
The intrinsic radio quiet nature of the AGN in Mrk 231 is not directly apparent because of optical reddening 
and Doppler boosting. 
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solid angle subtended by the BAL outflow. Mrk 231 is a LoBALQSO (a BALQSOs that 
exhibits low ionization resonant absorption troughs). Mrk 231 has broad absorption in a 
variety of low ionization species, Mg II A 2798, the well known Na I D doublet, He I A388 9 
I, Ca II K A3934 and Ca II H A3968 JSmith et all Il995l : Ihipari et allll994l : iMartel I Il997h . 
Thus, Mrk 231 displays the full panoply of quasar related phenomena, broad emission lines, 
powerful thermal luminosity, UV broad absorption lines and a radio jet. Therefore, as a 
consequence of its proximity to earth, Mrk 231 can provide an optimal opportunity to study 
the interaction of these physical elements by means of high resolution VLBI observations on 
sub-parsec scales. 

Earlier 15 GHz self- calibrated VLBA (Very Long Baseline Array jfl observations reported 



Ulvestad et al Jl999al jbl) with sub mas resolution revealed a double radio source separated 



m 

by approximately 1 pc. We initiated a program to look for more sub-structure and component 
motion at higher resolution by means of phase-referenced observations at 22 GHz and 43 GHz 
VLBA as discussed in section 2. No further resolution into sub-components was achieved 
even at 43 GHz. The two archival epochs showed that the nucleus was very dynamic and 
set a lower bound on the time var iability brightness tem perature of Tf, > 10^" K, below the 
inverse Compton hmit, ~ 10^^ K (lUlvestad et allll999a| ]bl). Thusly motivated, our program 



of observations attempted short time interval sampling (3 months between the epoch 2006.07 
and 2006.32) in order to elucidate this highly dynamic behavior. Mrk 231 did not disappoint, 
there was an enormous flare between these two epochs. Section 5 is dedicated to analyzing 
the rapidly changing core. We use our data and the adopted cosmological model to show 
that the strong flare will violate the constraints associated with the brightness temperature 
limit unless relativistic Doppler beaming is invoked. 

We also explore the steep spectrum secondary component of the nucleus of Mrk 231 in 
section 6. It is mildly variable and it is just less than a par sec from the core. Based on 
the nature of the relativistic core, it is natural to wonder if there is any relativistic motion 
of the secondary. In section 3, we determine that the secondary is either not moving at 
all or is moving at velocities far below relativistic values. Thus, the more relevant question 
becomes what is constraining the secondary from having such motion. Section 6 uses our 
VLBA observations and archival low frequency observations to model the apparent free-free 
absorption of the secondary emission. We argue that this same gas is conflning the secondary 
by ram pressure and that this gas is likely one and the same as the BAL wind. 



■^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



-4- 



The Observations 



High frequency (15, 22 and 43 GHz) observations were conducted at epoclis 2006.07 
(Project ID BP124A) and 2006.32 (Project ID BP124B) with all ten antennas of the VLBA. 
In addition, data from 2000.02 (Project ID BU013E) at frequencies of 8.4, 15 and 22 GHz 
and 8.3 GHz data from 2006.67 (Project ID BA080B) were retrieved from the VLBA public 
archive. The data were correlated on the VLBA correlator in Socorro, NM and calibrated 
with NRAO's Astronomical Image Processing System (AIPS) software packa ge. In addition 



to the standard calibr ation procedure for VLBI polarization observations (jCotton I Il993 



Leppanen et al. Ill995l ). an att empt was rnade t o correct for amplitude losses due to opacity 
using the AIPS task APCAL Leppanen I Il993h . 



At all epochs and frequ encies, except 2006.67, the data were initially phase-referenced 
( iBeasley and Conwayl 119951 ) to J1302+5748. At frequencies below 43 GHz the data were 
then phase self-calibrated, but at 43 GHz there was insufficient signal-to-noise to allow this. 
The phase reference source was imaged to provide the best possible model for the phase 
referencing procedure, but proved in any case to be a very close approximation to a point 
source. 

Figure [T] captures the images of the 8.4 GHz observations from 2000.02 and 2006.67. 
Figures [2] and [3] are the images at 15 and 22 GHz, respectively, for all three epochs and 
Figure m is the 43 GHz images from 2006.07 and 2006.32. AH images of Mrk231 presented in 
this paper are made with Briggs' robust weighting scheme (lBriggslll995l ). with the robustness 
parameter set to zero, with the exception of the 43 GHz image at epoch 2006.32 where natural 
weighting is used. All images are centered on the position of the southwestern, secondary, 
component at rectangular coordinate (0 mas, mas). 

Each of the images shows a pair of unresolved components separated by about 1 mas 
that vary in brightness, but do not vary significantly in separation. The NE component is 
the most variable and has the flattest spectrum and it is assumed to be associated with the 
AGN nucleus. 

Gaussians were fltted to the pos i tions of the two components using the DIFMAP soft- 
ware package (IShepherd et al. 1 1 19941 . Il995l ) and the resultant parameters are presented in 
Tables [T] and [21 Errors in the flux density in Table 1 arise from two sources: the absolute flux 
scale error, and an error due to uncert ainties in the model flt. The method used to determine 
the model fit part of the error is from iBritzen et all (120071) . The a bsolute fiux sc ale errors for 
the VLBA are estimated to be ~ 5% at 8 GHz Juivestad 1120081) and 15 GHz (IHoman et al 



20021), and ~ 7% at 22 GHz (IHoman et alll2002f ). At 43 GHz we estimate the error on the 
absolute fiux density scale to ~ —8%/ -|- 15% - the non-symmetric error arising from the fact 
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that some flux density is likely lost due to the inability to phase self-calibrate this source 
at that frequency. An analysis of the amplitude corrections derived from self-calibration of 
various cahbrator sources indicates that these observations were not affected by any unusual 
amplitude errors that would lead to larger errors than the values quoted above. 

Based on the inverted spectrum from epoch 2000.02, we anticipated dramatic behavior 
on the smallest of scales. Thus, we proposed for two epochs of high frequency observations 
that were closely spaced in time, 2006.07 and 2006.32. The observations were phase refer- 
enced to J1302+5748 at 15, 22 and 43 GHz. We also gathered polarization data. There is 
no evidence that the flux density recovered at any of our frequencies was strongly dependent 
on the weighting used, indicating that we are dealing with two compact components with 
little, if any, diffuse flux on the range of spatial scales probed by these observations. 

The inability to phase self-calibrate at 43 GHz opens the possibility that some flux 
density may have been lost due to phase errors. In lieu of the undetected secondary calibrator, 
J1306+5529, one can use the secondary component to the SW in MRK231 as a crude 
amplitude calibrator. Note the flux density of the secondary was the same between epochs 
2000.02 and 2006.32 at 22 GHz to within a few percent in Table 1 indicating that it is 
not strongly variable. Therefore, we treat this component as approximately steady at high 
frequency and it can potentially serve as a crude calibrator for the core flux density at 43 
GHz. The flux density of the secondary (our crude high frequency amplitude calibrator) 
measured at both epochs at 43 GHz agrees to within 1.2 mJy (0.6 sigma). This supports the 
suggestion that the phase referencing with respect to J1302-I-5748 recovered the majority of 
the flux at 43 GHz, and any amplitude loss due to phase errors is likely to be similar at the 
two epochs. 

We note that in epoch 2006.07 bad weather and technical problems allowed us to recover 
data from effectively just 5.25 of the 10 array antennas at 43 GHz, so these data have sparser 
u-v coverage and an increased RMS noise level. Our 43 GHz images in Figure H] show partial 
resolution of the core at epoch 2006.07 with an extension at PA = 135° from NE to SW. 
Unfortunately, we must consider this detection unreliable due to the poor beam shape and 
the low SNR. We stress that the data from all 10 antennas were recovered at 43 GHz at 
epoch 2006.32. This higher quality data is used extensively in section 5. 

In Figures 5 and 6, we plot the spectrum of the core and the secondar y, respectively, 



for all three epochs of observation as well as the two archival data sets in lUlvestad et al 
Jl999al lb[). 



The core spectra in Figure 5 are highly variable in magnitude. All the epochs reveal an 
optically thick core that is indicated by flat or inverted spectra in a finite band of frequency. 



Fig. 1.— VLBA 8 GHz images of Mrk 231 from epoch 2000.02 (left) and epoch 2006.67 
(right). The contours start at 0.7 mJy/beam and increase in factors of 2 to a maximum of 
89.6 mJy/beam. 
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Fig. 2.— VLBA 15 GHz images of Mrk 231 from epoch 2000.02 (left), epoch 2006.07 (middle) 
and epoch 2006.32 (right). The contours start at 1.8 mJy/beam and increase in factors of 2 
to a maximum of 57.6 mJy/beam. 
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Fig. 3.— VLBA 22 GHz images of Mrk 231 from epoch 2000.02 (left), epoch 2006.07 (middle) 
and epoch 2006.32 (right). The contours start at 1.8 mJy/beam and increase in factors of 2 
to a maximum of 57.6 mJy/beam. 
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Fig. 4.— VLBA 43 GHz images of Mrk231 from epoch 2006.07 (left) and epoch 2006.32 
(right). The contours start at 2.2 mJy/beam and increase in factors of 2 to a maximum of 
17.6 mJy/beam. 



Table 1 


: High Frequency VLBA Observations Data Summary 
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S.4 CIIx Flux 
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Density (nijy)/ A, T" 


Doiir,it\- (nijyj/ P'' 


Density (lii.7y)/Vi P'' 


Density (lii,]y)/% P" 


NE 


1996.93 


18 ± 2 


17 ±2 






NB 


1998.71 




44 ± 3 






NB 


2000.02 


29.1 ± 1.6/ < 1.0 


31.3 ± 1.8/ < 1.9 


26.7 ± 2.0 < 3.9 




NB 


2006.07 




43.1 ± 2.2/ < 1.9 


29.4 ± 2.2/ < 3.4 


ZO.ljl^-g/ < 13.7 


NB 


2006.32 




68.6 ± 3.6/ < 1.2 


74.8 ± 5.3/ < 0.7 


21.3+11/ < lO-l 


NB 


2006.67 


67.0 ± 3.9 /... 








SW 


1996.93 


114 ± 11 


51 ± 3 






sw 


1998.71 




60 ± 3 






SW 


2000.02 


126.2 ±6.3/ < 0.24 


71.2 ± 3.7 < 0.28 


36.5 ± 2.7 < 2.9 




sw 


2006.07 




69.6 ± 3.6/ < 1.2 


43.1 ± 3.1/ < 1.3 


S-Sjij I/ < 28.1 


sw 


2006.32 




66.8 ±3.4/ < 1.3 


40.0 ± 2.9/ < 9 


S-eto.g/ < 26.1 


sw 


2006.67 


127.6 ± 6.3 /... 







"The polarization is a 3 cr upper limit 
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Table 2: Component Separation 



Date 


8.4 GHz " 






15 GHz » 




o-2° 


22 GHz " 






43 GHz 






1996.93 








1.081 


0.03 


0.172 














1998.71 








1.162 


0.004 


0.174 














2000.02 


1.038 


0.004 


0.126 


1.173 


0.003 


0.139 


1.219 


0.005 


0.128 








2006.07 








1.009 


0.006 


0.172 


1.093 


0.005 


0.132 


1.222 


0.013 


0.142 


2006.32 








1.172 


0.002 


0.089 


1.193 


0.002 


0.059 


1.216 


0.010 


0.103 


2006.67 


1.124 


0.016 


0.267 





















"Measured separation at the designated frequency in mas 

^(Ti is the error computed using the technique of ICondon et all (|l998l ). see equation (1) 



^(72 is the error computed using the technique of Lobanov ( 20051 ). see equation (2) 
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Fig. 5. — The core spectra of Mrk 231 at six different epochs of VLBA observation. Notice 
that the data point for 1998.71 overlaps the 15 GHz point at epoch 2006.07. 
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Clearly the most powerful epoch is 2006.32 which we explore in detail in section 5. In contrast 
Figure 6 shows that the secondary is optically thin, steep spectrum at all epochs. The data 
are supportive of a secondary that is mildly variable, characterized by changes of < 10% 
(< la) on time scales of 0.25 years, ~ 15% (< 2cr) on time scales of 6 years and ~ 30% 
on time scales of 9.4 years. Notice that the spectra show significant curvature, reflecting a 
deficit of high frequency flux relative to a power law extrapolated from lower frequency. 



3. Component Motions and Positions 

An important part of this analysis is to try to determine the motion of the secondary 
relative to the core over the 9.4 year history of VLBA observations that were capable of 
resolving the compact nuclear double (i.e ., observations at > 15 GHz). The entries in Tables 
1 and 2 include the historical data from lUlvestad et all (jl999al Jbl) for this reason. It is clear 
that there is very little component motion. Therefore, a careful analysis of the errors that 
arise in the determination of the component separation is critical. The displacements between 
the core and secondary are entered in Table 2. Those entries in the first five rows have errors 
computed by a method that is based on the error in discerning the position of t he centroid of 
an elli ptical Gaussian fit to the data in the presence of background noise (e.g 
(Il998|)), 



Condon et al 



AO. 



Sep 



0.5 



'beam 



SNRr 



9 



beam 



-Y 

'S'A'^-Rsecondary / 



0.5 



where A^sep is the error in the estimate of the separation in mas, ^beam is the beamwidth in 
mas and SNR is the signal to noise ratio. Notice that the component separations in Table 



( Lobanov 


1998; 


Koniel 


1981) 



separation between epochs a common frequency should be used. The lone common frequency 
of observation for all epochs is 15 GHz. 

However, the method described by equation (1) does not account for core subcompo- 
nents that flare on scales below the resolution limit of the beam without any actual motion 
occurring. If a subcomponent flares, the centroid of the elliptical fit of the core can shift 
in position and give the illusion of secondary motion. This is a concern for Mrk 231 since 
Figure clearly indicates the presence of strong core flares that are not resolved by VLBI. For 
example, the core flux density flares at 1998.71 relative to the flux density of the core at 
1996.93 and the components appear to separate by 0.08 mas at 15 GHz according to Table 
2. This can be described by the base of the jet at the far eastern edge of the core being the 
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10 20 
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Fig. 6. — The spectra of the secondary component of the nucleus of Mrk 231 at six different 
epochs of VLBA observation. 
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source of the flare with no actual secondary motion. The position of the elliptical fit just 
shifts to the east by about 0.08 mas, well below the resolution of the beam at 15 GHz or 
22 GHz where the flare is detected. Thus, we need an estimate of the error in component 
separation that incorporates the maximum size subcomponent flare in the core that can not 
be resolved with VLBI (equivalently, the minimum size subc omponent that can be resolved 



by VLBI). Comparing the maps of the compact double in lUlvestad et all (jl999al ) at 15.3 
GHz (resolved), 8.4 GHz (partially resolved) and 5 GHz (unresolved), it appears that the 
uncertainty posed by unresolved subcomponents in the core will be about 1/2 - 1/3 of the 
beamwidth projected along the direction toward the secondary. Fortuitously for the com- 
pact double in Mrk 231 this direction is close to the minor axis of the beam. In general, 
the minimum size subcomponent that is resolvable by VLBI will depend on the SNR, so we 
look for a more general and quantitative formulation. To this end we invoke the technique of 



Lobanov I (120051 ) that considers the ability of an interferometer to resolve components mod- 
eled by elliptical fits in presence of noise assuming a rectangular power law interferometric 
visibility sampling function. 



22-/3/2 



lim 



TT 



Tiab In 2 In 



/ SNR \ 
\SNR-l J 



0.5 



(2) 



In equation (2), 6'iim is the minimum size subcomponent in mas that can be resolved by the 
interferometer. In our calculations we have the weighting of the data, /? = corresponding to 
unweighted data. The major axis and minor axis of the beam in mas are a and 6, respectively. 
The values of that we derive from equation (2) range from 6/4 - 6/2 for the observations 
in Table 2. We consider, an estimate of the uncertainty in the core position due to 
unresolved flaring subcomponents. The secondary does not flare, so we do not consider this 
a likely source of error in the position of the secondary. We add this error to the error 
introduced in equation (1) in quadrature (6'iim is the dominant contributor). The results 
computed by the formula, (2), are listed for the separation error on rows 6 - 10 of Table 2. 
They are large compared to the variations in the component separation between any pair of 
epochs. We conclude that the observations in Table 2 are not precise enough to determine 
the rate of component separation, except as an upper limit determined by On^. The VLBA 
observations indicate that the separation in components changed in 9.4 years (from 1996.93 
to 2006.32) by 

= 0.091 ±0.094 mas. (3) 

In arriving at the result given by equation (3), we noted that the core at epoch 1996.93 was 
at an extreme historic low. Thus, we do not consider core flares relevant in this quiescent 
state. Therefore, we use the error derived from equation (1) that appears in the first row of 
Table 2 as an accurate assessment of the error in the measured component separation at this 
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epoch. At 2006.32, there is clearly a flare, so we use the 15 GHz entry derived from equation 
(2) that appears in row 10 of Table 2 to estimate the error in the component separation at 
this epoch. The secondary separation in equation (3) equates to a time averaged secondary 
separation velocity, Vadv 

Vad^ = 0.026c ± 0.027c . (4) 

Our data is consistent with no apparent motion of the secondary relative to the core. It should 
be remembered that Vadv is the projection of the separation velocity of the components onto 
the sky plane. 



4. The Core Structure at 2006.07 

Table 2 indicates that the component separation increases with observing frequency at 
all epochs. This is expected in optically thick cores as the r = 1 surface for synchrotron 
self absorption moves outward relative to the true core position with decreasing frequency 



( iKonigl Ill98ll ). In order to ascribe a quantitative measure of this phenomena one must resort 



to highly model dependent estimates that involve m uch speculation about many unknown 



physical quantities that parameterize the putative jet fiLobanov Ill998l ). However, notice that 
this variation is particularly pronounced at 2006.07. Furthermore, from Figure 5, the core 
is actually optically thin between 15 GHz and 22 GHz, yet there is an 84 /ias differential in 
the component separations. In this section we study the ramiflcations of this to see if we can 
make some quantitative statements with a minimal amount of model dependent assumptions. 

We begin by overlaying the component positions on the 43 GHz radio image at epoch 
2006.07 in Figure 7. The main issue is how does one register the image centers at different 
observing frequencies. This was not designed as an astrometric experiment and we do not 
have a good handle on the positional jitter that is induced by fluctuations in the tropospheric 
delay. 

Ostensibly there is no obvious way to place errors on the component positions in Figure 
7. However, notice that the secondary components are fairly close in position. In fact we 
can argue that they should be even closer than this based on the steep spectrum in Figure 
6. We seek to register all of the epoch 2006.07 radio images relative to the 22 GHz position 
of the unresolved secondary at 2006.07. Let S22 be the region containing the plasma that 
emits the bulk of the 22 GHz secondary emission. For any reasonable ensemble of energetic 
particles, there must also be emission at 43 GHz and 15 GHz. First, consider the 43 GHz 
emission. The spectral index in Figure 6 between 22 GHz and 43 GHz is incredibly steep, 
all = 2.2, where the spectral index is deflned in terms of the flux density by the convention. 
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Fig. 7. — An overlay of the frequency dependent component positions on the 43 GHz image 
at epoch 2006.07. The errors in positions are naively computed using equation (1). Only 
the 43 GHz points have significant error bars. A more accurate error assessment method is 
shown in Figure 8 and discussed in the text. 
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Siy ~ z/~". Between 15 GHz and 22 GHz the spectral index is significantly fiatter, but still 
very steep, 0:22 = 1-2. Thus, unless there is an incredibly abrupt spectral break above 22 
GHz and a second weak 43 GHz component, virtually all of the 43 GHz radiation from the 
secondary would be emitted from the particles in 1^22 ■ Therefore, we consider the assumption 
that the 43 GHz emitting electrons are co-spatial with the 22 GHz emitting plasma is most 
likely valid. We can not make the same argument for the 15 GHz emitting electrons. Since 
the spectrum is so steep between 43 and 22 GHz and the spectrum above 22 GHz is also 
steep, we make the conservative assumption that the 15 GHz fiux density from S22 is at 
least as big as the 22 GHz fiux density, 43.1 mJy. With this conservative assessment, we can 
estimate an upper bound on the error in the 15 GHz secondary position if it is re-registered 
co-spatially with the 22 GHz secondary position. Namely, if the remaining 26.4 mJy of 
residual 15 GHz fiux resides primarily in a disjoint component, S15, this would create the 
maximum offset in the true absolute position of the 15 GHz component p osition relative 



to S22 . The 15 GHz secondary is unresolved in Figure 2, and according to lUlvestad et al 



(120051 ) and our own experience discussed in section 3, if S15 was half a beam width away 
from S22, the secondary would be partially resolved in the radio map at 15 GHz. There is 
no partially resolved component in Figure 2. Thus, we can set an upper bound on the error 
in registering the 15 GHz secondary position to the 22 GHz position, by positing that 26.4 
mJy of fiux density arise from an optically thin component, S15, that lies one half of a beam 
width from S22 (which accounts for the remaining 43.1 mJy of the 15 GHz fiux density of the 
secondary). Then taking the fiux density weighted average of the component positions (i.e., 
offset from S22 = (26.4/69.5)(6'/2)) yields a conservative upper bound on the offset between 
the true 15 GHz position and the 22 GHz position. The resultant upper bound on the error 
in the re-registered 15 GHz position is 0.11 mas and 0.018 mas in the x and y directions, 
respectively. Figure 8 is the Figure 7 overlay re-registered to the 22 GHz secondary position 
(i.e., all three of the secondary positions agree). The errors on the re-registered 22 and 43 
GHz core data points in Figure 8 arise from equation (1). The error on the 22 GHz core 
position is too small to be seen in this figure. 

The re-registered map in Figure 8 makes the component positions appear almost linear. 
We believe that our assignment of errors justifies further analysis of the frequency dependent 
component separation data. Apparently, we have uncovered sub-structure on scales smaller 
than the 43 GHz beamwidth. But, unlike the optically thick case in which the core position 
movement refiects the location of the frequency dependent r = 1 surface, the 2006.07 core 
spectrum is optically thin below 22 GHz and we are revealing direct data on the distribution 
of the synchrotron emitting gas. In particular, the spread in the 15 GHz and 22 GHz core 
positions in Figure 8, refiects the fact that there are significant amounts of 22 GHz emitting 
plasma to the east of the plasma that emits the 15 GHz data. In fact, the spectrum in Figure 
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Fig. 8. — An overlay of the frequency dependent component positions on the 43 GHz image 
at epoch 2006.07. This figure is similar to Figure 7, except that the components are re- 
registered relative to the 22 GHz secondary position. This tells us how the map center at 
one frequency is related to the map centers at the other frequencies. To get an absolute 
registration within the 43 GHz radio image, we place the 3 coincident secondary positions 
at the location of the best fit Gaussian to the 43 GHz secondary flux density. The A and B 
are labels that show the component locations within the theoretical two component model. 
Component A is synchrotron self absorbed and component B is very steep spectrum. 
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5 and core component positions are well represented by a two component model. Component 
A is located ~ 40 fias to the northeast of the 43 GHz core position, along a line connecting 
the three core position components and component B is ~ 40 /xas to the southwest of the 15 
GHz core position along the same axis. Component A is synchrotron self absorbed with a 
peak flux density of ~ 27 mJy at ~ 32 GHz. Component B is steep spectrum with a ~ 2.4. 
We only introduce the model to give a qualitative feel for the underlying substructure not 
to fine tune the parameters. This model explains the inflection point in the SED at epoch 
2006.07 inferred from Figure 5 and predicts ~ 3 mJy at 43 GHz in component B, similar to 
what we see in the southwest core extension in the radio image in Figure 4. Again, we caution 
that this feature might just be phase noise. Without a higher signal to noise detection, we 
can not say more. However, the spread in component positions and the SED are consistent 
with there being some emission in this region. 

The reason for pursuing this two component model in Figure 8 is not to make precise 
statements about the components that are not directly detected, but to explore the more 
general implications of a core substructure at 2006.07, composed of an optically thick com- 
ponent and a steep spectrum component, that is required to flt the SED. To this end, we 
plot the SED of the theoretical components A and B relative to the SEDs at 2006.07 and 
2006.32 in Figure 9. 

The first deduction is that the flare at 2006.32 can not be created by the simple adiabatic 
expansion of the fiaring core component A at 2006.07. If it were simple adiabatic expansion 
of a flare that is peaked at 32 GHz at 2006 .07, then the optically thin flux at 43 GHz should 



be much weaker at 2006.32, (jMoffetlll975l ). but it is actually stronger. Thus, the flare at 



2006.32 requires a large injection of energetic particles and energy after 2006.07. 

The second important implication is that component B is very steep spectrum above 15 
GHz. Figure 9 shows that there is some mechanism that causes the high energy spectrum 
to decay rapidly in the core of Mrk 231. All the components are extremely steep in the 
optically thin region, thus there is a very efficient cooling mechanism that preferentially 
cools the particles with the highest energy first. The most common mechanism in AGN that 
fit this prescription are inverse Compton cooling from a strong radiation field or synchrotron 
emission from a strong magnetic field. The cooling is so efficient that by the time that an 
ejection from the nucleus is far enough away so that it could be resolved by 43 GHz VLBA, 
it is generally too weak to be detected. Even with all 10 antennas working at peak efficiency, 
two good calibrators at 43 GHz and good weather, partial resolution of the ejecta might be 
the best one can do when studying fiare evolution in the nucleus of Mrk 231. 

The next section is the main topic of this study. We explore the formation of the 
powerful fiare at 2006.32. Furthermore, we try to get an understanding of the rapid cooling 
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THE CORE SED in the TWO COMPONENT MODEL 
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Fig. 9. — Comparison of SEDs of the theoretical two component model of the core at 2006.07 
with the SED at 2006.32 
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mechanism that causes drastic temporal changes in the emissivity of the ejected flares. 



5. The Core Flare at Epoch 2006.32 

The data in Table 1 indicates that the NE component is highly variable and is likely 
associated with the AGN nucleus. The two epochs of the observations, 2006.07 and 2006.32, 
reveal a strong flare at 22 GHz that increased the flux density by 45 mJy in just 91 days. The 
inverted spectrum from 15 GHz to 22 GHz juxtaposed to the steep spectrum, a = 1.8, from 
22 GHz to 43 GHz is extremely unusual for a radio source and provides a strong constraint 
on the nature of the flare that we exploit in this section. The steep spectrum above 22 GHz 
is inconsistent with free-free emission unless the temperature of the gas is < 10 ° K which is 
unlikely within a few pc of a quasar nucleus. Thus, we consider synchrotron radiation as a 
more likely source of the emission. 

There are two possible explanations for the inverted spectrum near 22 GHz. First, the 
spectral luminosity, Ljy, might be seen through a free-free absorbing (FFA) screen. For a 
powerlaw synchrotron spectrum seen through a FFA gas with number density of free elec- 
trons , A"fi, temperature, T, and thickness in parsecs, Lpc, we adopt the following expressions 



from IWalker et all fl2000[ ) , 



k{u) = 9.8 X 10-%cA^eT-i-V-2[i7.7 + ln{T^-^u-^)] . (6) 

We have only three data points in Table 1, the flux densities at 15 , 22 and 43 GHz of the 
core. If we rewrite the free-free opacity as k{FFA) = KoI^~^\11-1 + ln(T^'^z/~^)], we can 
separate out the frequency dependence, where the coefficient = 9.8 x lO^^LpcA^^fT"^-^. 
We assume that the emission region is a homogeneous sphere ejected along the jet direction 
toward the secondary. This is certainly more accurate than a continuous jet for emission 
that shows no evidence in Figure 5 or Table 1 of being quasi-steady. The lack of protracted 
periods of quasi-steady flux is in contradiction t o what would be expected from a large 



smoothly varying structure, like the jet models in iKonigl I (jl98ll ). Every indication is that 
the emission around 2006.32 is in the form of a large impulsive burst on the background 
of some low luminosity feature (perhaps the continuous jet). The episodic flaring nature 
of the core indicated in Section 4, Figure 5 and Table 1 is consistent with large ejections 
of highly energized plasma, with the epoch 2006.32 being the most energetic. The epoch 
2006.32 plasmoid is likely oblong and irregular with an inhomogeneous density, but without 
any prior knowledge of these geometric properties, we use a spherical solid as a zeroth order 
approximation. With this assumption of a homogenous sphere, the coefficient is a constant 
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in our calculations. We also assume for now that the covering factor, / = 1. For every value 
of T, we have three variables in equations (5) and (6), Kq, a and the normalization Lg. 
Thus, we can simultaneously solve for the three variables as three equations (the 3 measured 
flux densities converted into spectral luminosities in equation (5)) with 3 unknowns (the 3 
variables) and the resultant will fit the three flux density data points perfectly. We find the 
unrealistic result that a ~ 3.1 for the broad range of temperatures lO^K < T < 5 x lO^K. 
The extremely steep background spectrum is a consequence of the fact that the observed 
spectrum is very steep just beyond the spectral break, a ^ 1.8 from 22 GHz to 43 GHz. If it 
is the FFA absorbing gas that creates the spectral peak between 15 GHz and 22 GHz in Table 
1 then it necessarily has a strong effect on the spectral slope at 22 GHz. In particular, since 
k{FFA) is a function of u, the attenuation is highly significant at 22 GHz. Consequently, 
the background synchrotron spectrum at ^ 22 GHz would have to be appreciably steeper 
than indicated by the observed value of a = 1.8. 

Secondly, the spectral luminosity might be a consequence of radiative transfer through 
a synchrotron self absorbed (SSA) plasma. The strategy is to perform the final calculation 
in the plasma rest frame using known variables from observation. We designate the observed 
quantities with a subscript, "o", in the following expressions. Taking the standard result for 
the S SA attenuation coeffic i ent in the plasma rest frame and noting that u = Uo/S, we find 
from [Reynolds et al I (119961 ): iGinzburg and Syrovatskiil (Il969l ). 



9{n) 



r[(3n + 22)/12]r[(3n + 2)/12]r[(r2 + 6)/4] 



r[(r2 + 8)/4] 

The Doppler factor, 5, is given in terms of F, the Lorentz factor of the outflow; /5, the three 
velocity of the outflow and the an gle of propagation to the hne of sight, 6] 5 = 1/[F(1 — 
(3 cos 9)] (ILind and Blandfordlll985l ). This equation derives from an assumed powerlaw energy 
distribution for the relativistic electrons, iV(7) = iVr7~", where 7 is the thermal Lorentz 
factor and the ra dio spectral index a = (n — l)/ 2. The radiative transfer equation was 
trivially solved in IGinzburg and Syrovatskiil (119691 ) to yield the following parametric form 



for from the SSA source. 



(1 



-^J.{u)R 



) 



(9) 



where R is the radius of the spherical region in the rest frame of the plasma. We can 
make the same decomposition that we made for k,{FFA) for the SSA attenuation coefficient, 
jj{SSA) = yUo^'o"^'^^"'' • If we assume that the source is spherical and homogeneous then we 
have three variables in (9), fio, a and Lg. Since we have three equations (the 3 measured 



- 20 - 



flux densities converted into spectral luminosities in equation (9)) with 3 unknowns (the 3 
variables), we can fit the three flux density measurements perfectly. We find that a = 2.26 (or 
n = 5.52). Although this is incredibly steep for a radio source it is still much more reasonable 
than what we found for FFA. Thus, we consider SSA a much likelier interpretation of the 
peak in the synchrotron spectrum than FFA. 

Thusly motivated, we proceed to detail the spherical , homogeneous SSA model of the 
core. The synchrotron emissivity is given in iTuckerl (119751 ) as 



= 1.7 X 10~^\47rNr)a{n)B^^-^''\4: x 10^^)' 



a{n) 



2VV3)r(^)r(3^)r(n±^) 



8v/^(n + i)r (i^) 



a(5.52) = 0.102 . 



(10) 
(11) 



We can relate this to the observed flux density, 5'(z/o), in the opticall y thin region of the 
spectr um (43 GHz) using the relativistic transformation relations from lLind and Blandford 
fll985h . 



(12) 



where -Dl is the luminosity distance and is evaluated in the plasma rest frame at the 
observed frequency. A second equation arises from our SSA fit to the data above. 



/i(z/ = 19.5GHz)i?= 1 



(13) 



Since, we already determined a = 2.26 then the pair of equations (12) and (13) with the 
expansions of equations (7), (8), (10) and (6.7) indicate that we have two equations in four 
unknowns, R, B, 6 and A^r- We explore this solution space by setting R and 6 fixed then 
solving (12) and (13) for B and Nr, numerically. The results are displayed graphically in 
Figure 10. We solve for the actual number density, N, from the coefficient A^r by integrating 
the powerlaw distribution for F = 5.52 down to 7 = 1, A^ ^ Ar/4.22. We justify the use 
of 7 = 1 as a minimum of integration in our discussion of Table 3 and equation (18) below. 
We randomly chose four values of 6 that seem reasonable from VLBI studies of radio loud 
quasars, S = 1.00 ( non-relativistic) , 6 = 2.10 (mildly relativistic), ^ = 6.40 and 6 = 9.40 
(highly relativistic) (IKellermann et al II2004J : IVermeulen and Cohenlll994l ). For each value of 
S we found an exact powerlaw solution for both B and A^. The powerlaw for A^ is extremely 
steep. It is apparent that A^ and B take on extraordinary values outside of a narrow range 
of radii. These radii are smaller than the 43 GHz beam minor axis - roughly between about 
1/3 to 1/30 of the 43 GHz beam minor axis, so the emission region can not be resolved by 
VLBA. 
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The Number Density as a Function of Emission Region Size 
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The l\llagnetic Field Strength as a Function of the Size of the Emission Region 
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Fig. 10. — The number density (top) and the magnetic field strength (bottom) in a homo- 
geneous synchrotron self-absorbed spherical model of the 22 GHz fiare at epoch 2006.32 as 
a function of Doppler factor, 6 and spherical radius of the emission region, R. 
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Minimum Jet Kinetic Luminosity 
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Fig. 11. — The minimum kinetic energy flux, Q, for liomogeneous splierical source models 
cLS cL function of R and Doppler factor, 5. 
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It is straightforward to compute the total energy in a spherical volume from these 
solutions for B and A^, 



E{R) = J {UB + U,)dV=^nR^ / 



Kc2)(iVr7~"+i)d7 



(14) 



The justification for using 1 as the minimum value of 7 in (14) is given in our discussion of 
Table 3 below. The energy in the observers frame is Eo = 6E. From equation (14), we can 
derive a lower bound on the flare kinetic energy flux, Qmin{R) = SE{R)/tmax, where tmax = 
91 days is the maximum amount of time that it took the flare to grow (i.e., the separation in 
time between our two observations in 2006). We plot the resulting Qmin{R) values in Figure 
11 for each of the four values of 6. Each of the four plots of Qmin{R) has a minimum that 
we designate as the minimum kinetic energy flux for that value of delta, Qmin- To the left 
(right) of the minimum in each curve, the flow is thermally (magnetically) dominated. The 
minimum values follow an exact functional dependence on 6, logQmin = 44.36 — 1.81535. 
The properties of these four Qmin solutions are tabulated in Table 3. The first column in 
Table 3 is the value of 6 in the model. The next column is the minimum kinetic energy flux 
of the homogeneous spherical model of the flare. The next 3 columns are R, B and N for 
that Qmin solution. The synchrotron lifetime measur e d in years in the observers frame is 
calculated in column 6 using the formula from iTucken (119751 ) , 



t 



syn 



0.5 



sec 



(15) 



We compute this at 22 GHz in the observer's frame in order to see if the very steep spectrum 
from 22 GHz to 43 GHz in Table 1 can be explained by synchrotron cooling - spectral ageing. 
All of these lifetimes are short compared to 0.25 years between observations and since the 
flare probably expanded from smaller dimensions, B was probably significantly larger in 
early stages of the flare, so the synchrotron cooling rate was probably larger than indicated 
by column 6 in the early stages of the flare evolution. Thus, we conclude that synchrotron 
cooling is a viable explanation of the very steep spectrum, high frequency tail of the core 
spectrum in Figure 5. 

The next column tabulates the brightness temperatur e, Th, if the core flux c a me fr om 
the spherical volume of radius, R, using the definition from lScheuer and Williams! (119681 ) 



n 2kbZ/2 



(16) 



The values approach the inverse Compton limit that is given in iMarscher et al I (119791 ) as 
(see also equations (19) and (20), below) 



(17) 
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The large T;, values in Table 3 indicate that significant cooling of the relativistic electrons 
from inverse self-Co mpton losses are working in conso r t with syiichrotr on cooling to create 
the steep spectrum ( iKellermann fc Pauliny-Tothlll969l : lReadheadlll994l ). 



The last column in Table 3 is an estimate of electron energy 7 at the peak frequency. 
We know tha t the peak is le ss than 22 GHz, so these values are an upper bound. Using the 
formula form iTuckerl (119751 ). we compute the electron energy at 22 GHz to form an upper 
bound, 



O'peak 



z/o(peak)/5 



(3 X 1065)0-5 ' ' 

These are very low values for a relativistic jet and suggest that the energy spectrum likely 
goes down to 7 = 1 (i.e., there are electrons that radiate as low as 20 MHz, but the emission 
is not seen due to SSA). The low 7 values are related to the fact that all these minimum 
kinetic energy fiux solutions are very highly magnetized. 

The most revealing entry in Table 3 is the high value of T^. The r atio of inve r se Co mpton 
luminosity, Ljc, to synchrotron Luminosity, L synch was calculated in iReadheadI (jl994j ) as 



^ synch 

G{a , z) - 

fop/ (5 

3.5 GHz 



TbG(a, z 



|l/5 



10 



12.22 



1 



T\,G{a , z^ 



1 1/5 



+ 



fop/ 5 



1 + ^)1 GHz 



X 



10 



fhigh 



12.22 



(19) 
(20) 



1 + z)v, 



op 



(1+a) 



/3(«) 



"'(1 + ^)' 



where Uop = 19.5 GHz is the peak in the observed SSA spectrum, Uhigh = 43 GHz is the 
highest frequency in the power law in the observers frame (typically taken as 100 GHz in 
the rest frame of the plasma). Our choice of 43 GHz is conservative and Li^jL synch might 
be slightly underestimated by this assump tion. This expression a l so req uires the coefficient, 
/3(a), that is given by equation (15c) of IScheuer and Williams! (119681 ). For the minimum 
kinetic energy fiux solution with 5 = 1 in Table 3, the data in Table 1 and equations (19) - 
(20) imply that L irJLsmch ~ 2 x 10^. The X-ray observations from 0.2 keV - 50 keV from 
Braito et all (120041 ) that indicate a total intrinsic X-ray luminosity of Lx ~ lO^^ergs/s from 
the sum of accretion disk corona luminosity and the jet luminosity. Combining this with our 
measured radio luminosity of the core in Table 1, we observe that 



lO'^^ergs/s 



10^ 



(21) 



Lsynch lO^^ergs/s 
1 minimum kinetic energy fiux solution is inconsistent with the value of 



Thus, the 8 
Lie/ Lsynch determined by observation 
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One can investigate other possibilities in the 5 = 1 solution space to see if a solution 
that is compatible with equation (21) is viable. Essentially, R must be increased in order to 
lower Tfe. Increasing the spherical radius of the 6 = 1 solution so that Tf, satisfies equation 
(21) requires that R > 1.2 x 10^^ cm. According to Figure 11, this would correspond t o 
Q ~ 10^^ ergs/s, the power of the strongest known FR II radio jets (IWillott et al.lll999l ). 
These seems unlikely for a radio quiet quasar. 

Alternatively, if we look at the S = 9. 40 solution at R = 3.12 x 10^^ cm we get an 
intrinsic value of = Th/S = 3.7 x lO^'^ K ( jReadheadl 1 1994 ) . Inserting this into equations 
(19) and (20), implies that Lic/Lgynch ^ 10^ for this solution in agreement with the bound 
set by observation in equation (21). From Figure 11, this solution is magnetically dominated 
and Q ~ 3.5 x 10^^ ergs/s. This value of Q seems like a more reasonable value for a burst in 
a radio quiet quasars, on th e order of the long term time averaged Q in a strong FR I radio 
galaxy Jwillott et al.lll999h . 



In summary, we argue that the 5 
counts: 



1.00 solution is not physically reasonable on 2 



Tft is right below the Compton limit for the minimum kinetic energy flux solution. 
If this were true then there should be a strong intrinsic X-ray source (a nd it should 
diminish the radio emissivity as well) (IKellermann fc Pauliny-Tothlll969l ). Extensive 
deep observations and very har d X-ray observatio ns indicate that the intrinsic X-ray 
luminosity, Lx ~ 1 x lO^^ergs/s (IBraito et alll2004l ). This is actually at the low end for 
the intrinsic X-ray luminosity distribution from the corona of the accretion disk of a 
quasar which has the thermal luminos ity of Mrk 231 that was noted in the introduction, 
Lboi ^ X 10^^ ergs/s ( iLaor et all 119971 ). Thus, the strong intrinsic radio emission and 
the relatively weak intrinsic X-ray emission (by radio quiet quasar standards) is not 
consistent with near the Compton limit. Alternatively, R can be enlarged, but this 
presents different concerns as noted in the second point. 

The minimum flare kinetic energy flux is Qmin = 2.30 x 10*^^ergs/s. This is an Fanaroff 
Riley II (FR II) kinetic en ergy flux. The max imum FR I kinetic energy flux cutoffs at 
about Q = 5 X 10^^ ergs/s (IWillott et al.lll999l ). Also notice that the kinetic luminosity 
rises very quickly away from the minimum in Figure 11, so this is a very loose lower 
bound on Q. In order to alleviate the first concern on Tb, above, a significant departure 
from the minimum kinetic energy flux solution is in order as discussed in regards to 
equation (21) and therefore a marked increase in the flare kinetic energy flux. It seems 
unlikely (but not impossible) that a radio quiet quasar should have flares with a strong 
FR II kinetic energy flux. 
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Based on these two considerations, it seems difficult to reconcile the large kinetic energy 
flux with a Doppler factor of order unity in Mrk 231. Table 3 and the discussion of equations 
(19) - (21) indicates that a significant Doppler factor would ameliorate these concerns. It 
would also explain the gap between the core and the secondary in Figures 3 and 4. Small 
deviations in the line of sight can greatly alter the flux density since by equation (12), 
Sy ~ (5^'^^. Thus, small changes in the velocity or line of sight (jet wobble) can cause the 
jet emission between the core and secondary to be of such low surface brightness that it is 
too faint to be detected with VLBI. Similarly, a strong dependence on 5 is compatible with 
the extreme flares and spectral variations that appear in Table 1. It seems likely that the 
subparsec scale radio jet in Mrk 231 is highly Doppler enhanced, i.e. relativistic and viewed 
pole-on. This is not unexpected since numerous BALQSOs have been determined to be 
viewed in a pole-on orientation based on ra dio variability values above the Compton limit 
Jzhou et all [20061 : Ichosh and Punslvl boOTh . In particular, LoBALQSOs (those BALQSOs 
that exhibit low ionization resonant absorption troughs) seem inordinately common in these 
polar samples (1/3 - 1/2 of the polar BALQSOs). Mrk 231 is a LoBALQSO. 

There is strong corroboration in the literature for our deduction that Mrk 231 hosts a 
po lar relativistic jet. A 134ib38 mJy flux density variation in 1 day at 22.2 GHz was rep orted 



m 



McCutcheon and Gregory I (Il978l ). Using the methods of I Ghosh and Punslyl (l2007l ). the 
time variablity brightness temperature is = (12.4 ± 3.5) x 10^ ^ K. In order to avoid the 



inverse Compton catastrophe, the analysis in I Ghosh and Punslyl (120071 ) indicates that the 
line of sight to the jet is less than 9max = (25.6°) 



^+3.2° 
-2.2° 



We also note the fortuitous circumstance that 8.33 GHz VLBA observations were per- 
formed 0.35 years after the strong flare, epoch 2006.67 (BA80B). The historic high levels of 
8 GHz of core flux density at 2006.67 indicated in table 1 and figure 5 is supportive of the 
flare evolving to lower frequency 0.35 years later. Unfortunately we could not recover the 15 
GHz data at 2006.67 which would have told us much more. 



6. The Nature of the Secondary Component 



In this section we stu dy the interaction of th e secondary component with the surrounding 
medium. It was shown in lUlvestad et all (Il999al ). using low frequency (1.4 GHz - 15.3 GHz) 
VLBA observations that the compact radio double is seen through a FFA shroud. The FFA 
absorption creates a spectral turnover near 6 GHz. From Table 1, it is clear that the radio 
emission of the compact double is dominated by the secondary at low freque ncy. Thus, the 
emissio n measure (see equation (6), EM = f N^dl), that is calculated in lUlvestad et al 
(jl999al ) from this 6 GHz peak in the spectrum represents the shroud of gas that envelopes 
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the secondary and m ight not be indicative of the gas density surrounding the core. In 
Ulvestad et all (Il999al ). it was stated that EM ^ 1.2 x 10^cm~^pc produced the spectral 



turnover. However, it was not clear what temperature was assumed for the gas. Thus, we 
re-derive this result in a form that incorporates the temperature. 

In order to study the spectral turnover, we want to use the lowest frequencies possible 
in order define the shape of the spectral peak accurately. The compact parsec scale double 
discussed thus far in this paper appears as an unresolved "core" at low frequency. The 
"core" is part of a triple radio source at low frequency. Components to the north and 
south of the core are detected at 8.4 GHz and below. At 2.3 GHz, the resolution is not 
sufficient to avoid confusion between the "core" flux and flux from the north/south extended 
components. Thus, as a compromise, we chose data at a resolution matched to the 5 GHz 
observation. We flt the data using equation (5) and a background spectral index of a = 1.7 
for the synchrotron source. Our best flt used a covering factor, / = 0.98. The correspo nding 



solution for the emission measure is close to the value stated in lUlvestad et all (Il999al ). 

EM ^ 9 X lO^cm-'^pc (^j^r^ , 2 x 10^ K < T < 5 x 10^ K . (22) 



Considering that the secondary is surrounded by a dense cloud of free electrons as 
indicated in equation (22) and the secondary does not move appreciably according to Table 
2, even though the core jet is relativistic, it is natural to consider the possibility that the 
secondary is being braked to a slow speed by the surrounding gas that supports iVg. In 
physical language, is the pressure of the secondary in dynamic equilibrium with the ram 
pressure of the surrounding gas? To explore this possibility we flrst try to estimate the 
pressure within the secondary. 

To estimate the pressure, we flrst flnd the Q associated with the secondary and use the 
physical size to convert it to a pressure. Since the secondary is very steep spectrum at all 
epochs and is mildly variable, we consider it to be a compact analog of the radio lobes that 
are seen in classical double radio sources 0. Furthermore, this compact lobe is seen through 
the effects of a FFA shroud and we can not detect its low frequency emission, even though 



"^One might also think of this feature as a strong knot in a jet that disappears, then reappears as a diffuse 
low frequency l arge sc ale elongated structure (that is depicted by the low frequency VLBA observations in 



Ulvestad et all (|l999al) ) that is barely detected with 8.4 GHz VLBA at a PA 65° offset from the sub-pc scale 
jet direction. The very steep spectral index in figure 6, above 8.4 GHz (which equates to a background 
spectral index of a = 1.7 in the FFA fit above), seems more indicative of lobe emission than knot emission. 
For such a steep spectrum structure, the distinction of knot or lobe is probably irrelevant with respect to 
the application of estimate of the jet kinetic energy flux in equation (23) 
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it is plausibly > 1 Jy considering the 8.4 GHz flux density in Table 1. Perhaps th e most 



sophisticated estimate of jet kinetic energy flux was developed by lWillott et al.l (Il999l ) which 
incorporates deviations from the overly simplifled minimum energy estimates. The method 
also considers geometric effects, flUing factors, protonic contributions and low frequency 
cutoff. The long term t ime average o f the i et kinetic energy flux of a steep spectrum radio 



lobe, Q, is estimated in lWillott et al.l (119991 ): 



Q fa 3.8 X lO^^L^/^ergs/s , (23) 

where L151 is the total optically thin luminosity from the lobes in units of lO^^ergs/Hz/sr at 
151 MHz. The obstacle to implementing equation (23) is that it depends on the 151 MHz 
flux density. This is problematic since the 151 MHz flux from the secondary is blocked from 
our view by the FFA shroud. The lowest frequency measurement that we have in Table 1 is 
8.4 GHz and we need to extrapolate this value to 151 MHz. We know that the spectrum is 
very steep between 8.4 GHz and 15 GHz, a = 1.38, even with some diminished flux at 8.4 
GHz from FFA. There are very few (if any) known radio lobes that are this steep between 
151 MHz and 8.4 GHz. Thus, there is no justiflcation for extrapolating the 8.4 GHz flux 
density with a = 1.38. A more typical value for st eep lobe emission seems to be a ~ 0.8 



(IKellermann et al lll969l : lHerbig and Readheadlll992l ). Under this assumption, we used Table 



1 and equation (23) to estimate the long term time average kinetic energy flux of the putative 
jet connecting the core to the secondary, Q ^ 9.9 x lO^^erg/s. 

Even though this estimate of Q is crude, we can use it to estimate the pressure in the 
secondary. We can approximate the energy flux through the jet as Q ~ {l/4:)7iD^VadvU, 
where D is the diameter of the lobe, Vadv is the velocity that the lobe is advancing into the 
external medium and U = Ue + Ub is the total internal energy density. Since in a relativistic 
gas the pressure P/obe = (1/3)^^, we estimate 

40 

Plobe ^ „ . (24) 

If the lobe is conflned by the ram pressure of an external medium moving with a velocity, 
Vext, and density, rripNext, then we compute the pressure balance at the leading edge of the 
lobe to flnd 

Plobe ~ „ ~ mpNextiVadv " Vextf ■ (25) 

If we implement a flducial distance of Ipc and assume that the medium is fully ionized 
then by equation (22), Next ~ lO^cm"'^. We also get a value of D = 7.70 x lO^^cm from 
the angular size of the emission region, 0.31 mas (this was determined from model flts to 
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the 15 GHz VLBA data by lUlvestad et all (Il999bl ) and was verified by our fits in various 



epochs as well). We can then solve (25) with these parameters derived from the FFA fit and 
the VLBA data to get a qualitative understanding of the ram pressure confinement of the 
lobe. We explore equation (25) in two limits, first the enveloping gas is stagnant v^xt = 
which yields Vadv = 0.012c = 3490 km/sec. Secondly, we note that the distance from the 
secondary to the core is typical for the location of the LoBAL gas. The velocity of the 
Mgll absorption trough relativ e to the Ha emiss ion (which we take as an approximation 



to the quasar rest frame) from ISmith et all (Il995t ) can be equated with the velocity of the 



enveloping medium, Vext ~ 0.0173c, from which equation (25) implies that Vadv = 0.025c = 
7570 km/sec. These values are consistent with the slow advance speed indicated in Table 
2. The dynamic estimates correspond to a separation increase between the core and the 
secondary in 9.4 years of 0.041 sin ^ mas and 0.089 sin 6' mas, respectively (where 6 is the 
angle between the line of sight and the jet propagation axis). These values are consistent 
with our observations which by equation (4) bounds Vadv < (0.053/ sin 6')c, at the one sigma 
level. 

In summary, the relativistic jet is over-pressurized relative to the gas density of the 
FFA absorbing gas, so it advances. However, the FFA gas is very dense and the kinetic 
energy flux is modest, thus the jet can drill through the dense medium relatively slowly at 
a non-relativistic rate. Independent of the assumptions involved in deriving and evaluating 
equation (25), this qualitative conclusion is robust and is consistent with the observed VLBA 
component separation data in Table 2. 

This analysis raises the question, what is the medium that is confining the relativistic jet 
and braking its velocity to a virtual halt? Is it the putative dusty molecular torus, Vext = 0, 
or is it the BAL gas, Vext = vbal^ We argue that the FFA absorbing screen is more likely 
the BAL gas than the molecular torus based on the following points. 



The low frequency VLBI observations of iLonsdale et all (120031 ) led them to conclude 



that "the lack of 21 cm absorption or maser amplification against the compact core 
(which is the compact double with our higher resolution) shows that the molecular 
disk does not cover the central source." 

The high Doppler factor inferred from our modeling of the core flare indicates that 
the jet is viewed in a nearly pole-on orientation. Furthermore, the existence of BAL 
troughs in the spectrum of Mrk 231 indicates that our line of sight is through the BAL 
gas. Thus, the jet direction seems pass through the BAL gas on pc scales. 



In summary, the observational data is consistent with the jet ramming through the BAL 
gas and inconsistent with the jet plowing through the molecular disk. Thus, we proceed to 
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consider the BAL e; clS clS the FFA medium. 



If the jet is ramming through the BAL gas then the jet can be used as a remote probe 
of properties of the BAL region. Thus, MRK 231 might be uniquely suited for studying 
the BAL gas. We find that the properties of the gas derived in this section are consistent 
with pubhshed interpretations of BALQSO physics. The main piece of information is the 
emission measure in equation (22), that does not have a unique decompositi on in terms of 
the num ber density distribution iVe(r) and the extent of the FFA region, L. In lUlvestad et al 
(jl999al ). they chose L = Ipc and a constant A'^e = 1.1 x lO^cm"'^. We i nterpret this va lue 
of L within the context o f the polar models of BAL outflows flows in iPunsly Jl999al lb[): 
Proga and KallmanI (120041 ). Consider a conical outflow that experiences very little accelera- 
tion as might be the case in the LoBAL region (near UV lines) s ince most of the lin e driving 
acceleration is provided by far UV and EUV absorption lines (IMurray et allll995l ). In this 
scenario, by conservation of mass, Ne{r) = No/r"^ and we can integrate the emission measure 
in general form as 



EM 



(26) 



where the observed Lobs is the projection of L in the sky plane and 6 is the ang l e to th e 
line of sight to the radio jet axis (symmetry axis of the BAL outflow). In iPunslyl (Il999bl ). 
it was calculated that 6 = 15° was typical for a LoBAL or from Table 2, L = 3.75 pc and 
from equations (22) and (25), N(,{r = L) 4.9 x lO^cm"^ ( io^k )" "^^^ assume that the 

absorbing colum n is almost completely ionized {N^ ~ A^) in the BAL region, as deduced in 
HamannI (jl998l ). then we can estimate the total hydrogen column. 



N, 



H 



N{r)dr 5.6 x 10 



22 



T 



104 °K 



0.75 



cm 



(27) 



The total hydrogen col umn is sufficient to make the BAL absorption lines optically thick 
to resonant absorption (IMurray et all Il995l : iHamannl Il998l ) . Superficially, one would think 
that the value of Nh in equation (27) is incor npatible with th e BAL gas in Mrk 231 . This is 
actually not the case as discussed at length in iHamannl (119981 ): iHamann et al I (120031 ). It was 
found that many BALQSOs show absorption in PV AA1118, 1128 with an optical depth > 0.2. 
For such a rare species this is impossible, unless the CIV AA1548, 155 troughs are s e verely 
saturated. After studying many high resolution spectra of BALQSOs, iHamann et al I (120031 ) 
concluded that over half of the BALQSOs with deep far UV spectra showed significant 
PV AA1118, 1128 absorption. The implication was that these BALQSOs must have CIV 
A A 1548, 1550 troughs that are completely saturated with the bottoms of the troughs filled 
in by scattered light and un-attenuated (or mildly attenuated) lines of sight that are a 
consequence of partial covering of the central AGN by the predominant thick absorber. 
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After, numerically modeling the absorption region, iHamannl (119981 ) found that the CIV 
AA1548 1550 troughs are saturated as a consequence of absorption presented by a column 
that contained almost completely ionized hydrogen, lO^^cm"^ < < lO^^cm"^. The upper 
limit is set by the fact that the gas must be optically thin to Compton scattering. It was 
concluded that these large c olumn densities are lik ely common in BAL regions. Based on the 
analysis of iHamanru (119981 ): iHamann et al I (120031 ). the values of Nh determined by equation 
(27) are typical of many BALQSO absorbing columns. Furthermore, we remind the reader 
th at there is very little neutral hydrogen absorption towards the nuclear secondary as noted 



Lonsdale et all (120031 ). Yet by equations (22) and (27) there are large values of EM and 



m 

Nh which implies that the intervening hydrogen column density must b e almost c omple tely 
ionized as expected from the ionization r aodels of dense BAL winds in iHamanru (119981 ). In 
particular, the photo-ionization models of lHamannl (119981 ) yield a range of allowed ionization 
states: for a total hydrogen column density of = 6 x lO^^cm"^, the neutral hydrogen 
column density is quite small, 5 x lO^^cm"^ < Nhi < 5 x lO^^cm"^. 

From our analysis above, we have enough parameters to estimate the kinetic energy flux 
of the the BAL wind, Q{BAL) in terms of the parameters of the enveloping gas surrounding 
the leading edge of the secondary. First, assume that the BAL wind has a bulk velocity in 
the range of velocity spread given by the Mgll BAL absorbing gas and Nal D absorbing gas 
relative to the quasar rest frame, 0.0151c < vbal < 0.0173c. From the numerical models in 
Punslyi ( 1l999bl ). a polar LoBAL wind will have an opening angle of ~ 20° — 30°. Combining 
this opening angle and the BAL wind velocity with the 3.75 pc displacement from the nucleus 
and Ne ~ 4.9 x lO^cm"^ yields 2 x lO^^ergs/s < Q{BAL) < 6 x lO'^^ergs/s. It is curious that 
this is similar to the kinetic energy flux of the relativistic flare that was discussed in section 5. 
One can compare th is parametric analysis, with the theoretical polar BAL models in Table 
1 of iPunslyl (jl999bl ). First of all, for the numerical models, the kinetic energy flux in the 
BAL wind solutions is 7.2 x 10*^^ergs/s < QtheoreticniiBAL) < 2.3 x lO^^ergs/s, in agreement 
with what we deduced for MRK 231. However, the mass density in the BAL wind of Mrk 
231 is estimated to be a factor of 10 - 50 times larger than that which is predicted by mass 
conservation in the numerical models. A possible explanation of the additional gas density 
is that the polar BAL entrains gas as it propagates through the surrounding medium. This 
seems reasonable on two counts. First of all, there is the dusty, enveloping gas that partially 
obscures the nuclear region and reddens the spectrum. This gas distribution could be a 
consequence of the merger that is in progress with a nearby smaller galaxy to the south 
(lArmus et all 1 19941 ). So there is a potential for significant gas entrainment. Secondly, the 
BAL absorption is deep, but not particularly blue-shifted for a BALQSO, only about 5000 
km/s: this is consistent with a resistive drag resulting from mass loading of the wind. 
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Discussion 



In this article, we used liigli frequency VLBA observations to determine that the central 
engine of the nearby BALQSO Mrk 231, ejects relativistic plasma along a trajectory close to 
the line of sight. Furthermore, it appears that this jet is interacting strongly with the dense 
BAL gas at a de-projected distance ~ 3 - 4 pc from the radio core. 

The most striking finding was the strong 22 GHz flare that emerged from the core 
between epochs 2006.07 and 2006.32. The core spectrum, during the flare, is characterized 
by extreme gradients, namely it transitions from being very steep spectrum, a ^ 2 above 
22 GHz to being inverted a < at ~ 15 GHz. The magnitude of the flare combined with 
the three month interval between observations sets a lower bound on the time variability 
brightness temperature of T;, > 2 x lO^*' K, below the inverse Compton limit, ~ 10^^ K. 
However, all attempts in section 5 to model the high frequency peak of the spectral turnover, 
19.5 GHz, in combination with the steep spectral index above 22 GHz (a ~ 2) indicate that 
the flare is synchrotron self-absorbed and Tj, ^ 10^^ K, unless the flux density is Doppler 
boosted. Our Doppler boosted models indicate a kinetic energy flux, Q ~ 3 x 10^^ ergs/sec 
and an intrinsic (rest frame of the plasma) brightness temperature > 10^^ K. 

Our discovery of a relativistic je t in a radio quiet quasar does not stand alo ne. There 
are the aforeme ntioned BALQSQs inlGhosh and Punslvl (120071): IZhou et all (120061) as well as 
PG 1407+263 in lBlundell eTH J2003h and HI Zw 2 in lsrunthaler et all J200oh . Bv definition. 



most radio quiet quasars are weak radio sources, so there are not many observations with 
VLB! Consequently, we can not say if the relativistic jets that have been detected are more 
the exception or the rule. 

The compact nature of the jet emission in MRK231 and the strong interaction with 
the surround ing medium has many strong similarities to GHz peaked (GPS) radio sources 
(lO'DealEoOsh . In the interval between 2006.07 and 2006.32, Mrk 231 seemed to transform 
from a weak version of a GPS source (as it was in 1996.93) to a weak version of a rare 
class of radio sources kriown as high frequency peakers, but on a scale 1000 times smaller 
( lOrienti and Dallacasa I l2008l ). It i s weakly pol arized like GPS sources, but unlike GPS 
sources Mrk 231 is highly variable (lO'Deal Il998l ). We were unable to reliably detect any 
morphological changes associated with this change in character. It would be interesting to see 
if long term VLBA monitoring at 43 GHz and 22 GHz would reveal ejecta emanating from the 
core ~ a month after a flare. Another interesting aspect of high frequency VLBA monitoring 
would be to understand the flare growth and decay time scales. The information gleaned 
from this data would further constrain the physical model of section 5. This monitoring 
would undoubtedly reveal large time variability brightness temperatures and might confirm 
the presence of synchrotron cooling, i.e. is there a short lived 43 GHz precursor to a 22 GHz 
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flare. If tlie basic SSA model of the flare that was presented in section 5 is confirmed by 
these observations, a natural question arises: what is the origin of the strong magnetic field 
, ~ a few G, on scales approaching 10^^ cm? If the field propagates from the base of the jet, 
does it tell us something about the turbulent plasma inside the accretion disk? 

Another important finding was that the core was unpolarized even during a flare. This 
implies that the magnetic field is not ordered or there is a huge amount of weakly magne- 
tized gas around the central engine which depolarizes the synchrotron emission by Faraday 
rotation. The arguments in section 5, regarding seemed to favor magnetically dominated 
models. It is hard to understand how a strong field would not be ordered in a magnetically 
dominated plasma. This would seem to indicate a dense magnetized envelope surrounding 
the radio core. Perhaps this is the source of the seed field that is ultimately responsible for 
the strong magnetic field in the core flares. 

It is intriguing that in spite of the highly relativistically beamed Doppler core that the 
steep spectrum secondary (a > 1) is mildly variable and has very little apparent motion. 
This peculiar situation was studied in section 6. A standard estimate of the long term time 
averaged kinetic energy flux of the secondary reveals Q ^ lO'^^ergs/sec which is consistent 
with the core having episodic flares of Q ^ lO^^ergs/sec, if the two components are part 
of the same jetted beam of plasma. Thus, we established consistency between the analysis 
of the core in section 5 and the kinetic energy flux of the secondary that was estimated in 
section 6. We also noted that matched (to 5 GHz) resolution VLB A observations indicates 
that the secondary is seen through a shroud of free-free absorbing gas with an emission 
measure of ~ 10^cm~^pc. We argue that the steep spectrum secondary seems to be a radio 
lobe associated with the jet advancing into a dense medium (the jet is confined by ram 
pressure) that is also the source of the free-free absorption. The properties of the dense gas 
are consistent with a temperature, 10, 000 K — 20, 000 K, a displacement from the nucleus, 3- 
4 pc, and a total hydrogen column de nsity, Nh = 10^^ — lO^ ^cm"^. These gas parameters are 



often associated with the BAL wind (IHamann et al II2003I ). The implied kinetic luminosity 



of this BAL wi nd is wi t hin the narrow range of values that are predicted by the polar BAL 



wind theory in iPunslyi (jl999bl ) . The density is at least an order of magnitude larger than 
what is expected from such a wind, indicating that there is substantial gas entrainment as 
the BAL wind propagates outward from the nucleus. 

The jet in Mrk 231 might offer a unique opportunity to study the central dynamics of a 
radio quiet quasar. Not only is the quasar near to earth, but the jet is the brightest at high 
frequency of any radio quiet quasar jet. Thus, no other radio quiet quasar central engine can 
be explored with such high resolution (VLBI). Furthermore, now we know the orientation of 
the jet. The high resolution of 43 GHz VLBA allows us to peer into the inner regions of this 
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quasar, using the relativistic jet as a probe of the X-ray emitting gas, the X-ray absorbing 
gas, and the BAL wind. In summary, there are many potential discoveries that can be made 
with future VLBA monitoring coordinated with observations at other frequency bands. Our 
observations lead us to ask the the following questions 

1. Is there more evidence of superluminal motion, like apparent velocities of radio emission 
larger than the speed of light that can be detected? 

2. What is the maximum kinetic energy flux of the high frequency flares? 

3. Are these types of flares present in other radio quiet quasars, but with varying magni- 
tude? 

4. Why arc the flares so powerful, yet cool off so fast that MRK 231 never becomes a 
strong FR 1 radio source? More speciflcally, why do all the components in the the core 
and secondary seem to show such steep high frequency spectra? What is the powerful 
cooling mechanism? 

5. Is this cooling mechanism found in all radio quiet quasars and why does it not occur 
in radio loud quasars? 

6. Is the gas that entrains the secondary one and the same as the BAL gas? 

7. How fast do strong flares grow and what are the details of their rapid decay? 

8. Do the detected (but absorbed) X-rays come primarily from the same plasma that 
is responsible for the 22 GHz flares (i.e., the accretion disk X-rays are completely 
attenuated)? Do we see correlated flares in the X-ray and microwave bands? 

9. Since the microwave flare emission is from such small scales, it is nearly coincident with 
the scales of X-ray emission and absorption. Thus, can we expect a strong interplay 
between microwave flares and the X-ray spectral properties of MRK231? 

10. If the X-rays in Mrk 231 are predominantly inverse Compton emission from the mi- 
crowave jet, does the dense X-ray absorption column constrain the physical parameters 
of the enveloping medium that gets entrained by the propagating jet and thwarts it 
propagation to large scales? 

11. Does a strong microwave flare blast holes in the X-ray absorbing column allowing us 
to get a brief mildly attenuated glimpse of the X-ray emission from the accretion disk 
corona? 
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12. Do flares in the microwave correspond to a decrease in the X-ray emission of the 
accretion disk corona as has been hypothesiz ed as a universal phenomenon in black 



hole accretion systems (the "low-hard" state) (IMaccarone et al I l2003l : iMarscher et al 
2002h ? 



Based on these considerations that were raised by our observations, it seems clear that a 
coordinated observing program of VLBA and X-ray observations would shed light on many 
of these issues, in particular points 8 to 12. Generally speaking, coincident X-ray and VLBA 
monitoring is potentially a powerful method of studying any "disk-jet connection." 
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Table 3: Details of the Minimum kinetic energy flux Solutions 



s 


Qmin 


R 


B 


N 


^syn 




7 




ergs/s 


cm 


Gauss 


cm~^ 


Vo = 22 GHz 




peak 












yrs 


b 


c 


1.00 


2.29 X 10^4 


6.50 X 10^^ 


4.43 


3.35 X 10^ 


1.15 X 10-2 


0.86 


< 40.7 


2.10 


5.94 X lO^-'' 


4.75 X 10^^ 


2.69 


8.80 X 10^ 


1.68 X 10-2 


0.66 


< 36.1 


6.40 


7.86 X lO'*^ 


2.90 X 10^^ 


1.14 


1.81 X 10^ 


3.48 X 10-2 


0.47 


< 31.7 


9.40 


3.94 X 10^2 


2.47 X lO^'' 


0.86 


9.75 X 10^ 


4.38 X 10-2 


0.41 


< 30.1 



"measured in observers frame, see equation (15) 

''See equations (16) and (17). Based on 68.6 mJy at 15 GHz. 

'^see equation (18), fc'o(peak) < 22 GHz 



